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^- Abstract 

^ : 

<^ ! Prospects for future measurements of Deeply Virtual Compton Scattering at HERMES are 

00 ', studied using different simple models for parameterizations of generalized parton distributions 

I (GPDs). Measurements of the lepton charge and lepton beam helicity asymmetry will yield 

■ important input for theoretical models towards the future extraction of GPDs. 

» ■ 

O ; 1 Introduction 

oo , 

O . 

y—( . The study of hard exclusive processes in the Bjorken limit is now considered as a promising tool 
^ ] to gain new insight into details of the nucleon structure that cannot be studied with inclusive 
deep inelastic scattering (DIS). A unified theoretical description of hard exclusive and inclusive 
processes has been obtained through the formalism of generalized parton distributions JI], 0, |^ 
(see recent reviews in Ref.s P!B[1) which are also called skewed, off- forward or non-forward 



Oh: 

I 



^ . parton distributions. 

' An ordinary parton distribution represents the probability to find a parton with a specified 
^ ■ longitudinal momentum fraction x in the fast moving hadron and thus is summing over all 

■ partonic configurations containing such a parton. In contrast, GPDs represent the interference 
of two different wave functions, one with a parton having momentum fraction x + ^ and another 
one with a parton having momentum fraction x — ^. GPDs, besides the longitudinal momen- 
tum fraction variables x and C, (called skewedness parameter), depend on a third independent 
variable, the momentum transfer = (p — p'Y between initial and final nucleon states with 
momenta p and p', respectively. 

There are four different types of quark GPDs contributing to the simplest hard exclusive process: 
Deeply Virtual Compton Scattering (DVCS), ep — > ep'j. In the unpolarized distributions, 
H'^{x, ^, A^) and E'^{x, ^, A^), the quark helicities are summed over. The polarized distributions, 
H'^{x, ^, A^) and E'^{x, ^, A^), are responsible for the differences between right- and left-handed 
quarks. 

The generalized parton distributions combine the characters of both the ordinary parton dis- 
tributions and of nucleon form factors. On the one hand, in the limit A^ ^ 0, ^ — holds 

H'i{x,0,0) =q{x), ^^(x,0,0) = Ag(x), (1) 

where q{x) and Aq{x) are the ordinary quark number density and quark helicity distributions. 
On the other hand, the first moment of GPDs must satisfy the following sum rules. 
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j\xH''{x,^,A') = g\{A') 



1 



cixE^(x,e,A^) = h\{A'), (2) 



where (A^) and F|(A^) are the Dirac and Pauh form factors and g\{A^) and /i^(A^) are 
the axial- vector and pseudo-scalar form factors, respectively. 

In the above formulae the variable x is defined in the range (—1, +1) and negative values of it 
correspond to anti-quark distributions in the following manner: 

q{—x) = —q{x), Aq{—x) = Aq{x). (3) 

Two different regions exist for GPDs with respect to the variables x and ^. For |x| > ^ the 
GPDs are the generalizations of the ordinary parton distributions, while for |x| < ^ the GPDs 
behave like meson distribution amplitudes. 

The recent strong interest in GPDs was stimulated by the finding of Ji that the second 
moment of the unpolarized GPDs at A^ = is relevant to the spin structure of the nucleon 
since it determines the total quark angular momentum: 

J-1 r^dxx \h'^{x,C,A'^ = 0) + E''{x,C,A'^ = 0)\ . (4) 

The total quark angular momentum Jq decomposes as 

Jg = IaS + L, , (5) 

where AS/2 and Lg denote quark spin and orbital angular momentum, respectively. As AS is 
measured through polarized DIS experiments, a measurement of Jg through the sum rule Eq. 
(Q) in terms of GPDs provides a model-independent way to determine the contribution of the 
quark orbital momentum to the nucleon spin. Eventually even the contribution of the total 
gluon angular momentum Jg may become accessible through 

1 = Jg + J9- (6) 

First measurements of the DVCS lepton helicity asymmetry have been accomplished recently 
by HERMES § at 27.5 GeV and by CLAS || at 4.25 GeV. Several plans exist for further 
measurements of DVCS and of other hard exclusive reactions to accomplish a first insight into 
GPDs. At HERMES, an upgrade of the spectrometer with a recoil detector [0 will significantly 
improve the separation of exclusive events. The DVCS process has also been observed in 
e~^p collider experiments at DESY by ZEUS |^ and HI 0]. The DVCS cross section was 
measured and compared to the QCD-based predictions. Measurements of the lepton beam 
helicity asymmetry by HI and ZEUS will become possible in the near future when longitudinally 
polarized leptons will be made available also to the collider experiments at HERA. 
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The main aim of this paper is the evaluation of the anticipated statistical accuracy for future 
measurements of DVCS asymmetries at HERMES upgraded with the recoil detector. In section 
2 of the paper different versions of GPD parameterizations are discussed. The third section 
deals with the assessment of the expected size of the asymmetries, based upon different GPD 
parameterizations, and with the evaluation of their projected statistical accuracy. Finally, the 
conclusions of the paper are presented. 



2 Parameterization of Generalized Parton Distributions 

Two examples of GPD calculations are presently known in the literature. While bag model 
calculations [|I^ show a weak dependence of the distributions on the skewedness parameter, 
chiral quark soliton model calculations [jlT|, in contrast, show a strong dependence on ^. The 



common approach at the moment is to use a guess that satisfies general constraints on GPDs 



known from theory. This paper basically follows the ansatz proposed in Ref.s ||12| , |T^. Here, 
the dependence of the GPDs on is taken in a factorized form with respect to the other 
variables whereby satisfying Eq.(|^). Any scale dependence of the GPDs is neglected. 

In the simplest approach GPDs can be assumed to be independent of the skewedness parameter 
^. In the following, only u- and d-quark GPDs are considered to be non-zero. The function 
H, for example, is written as a product of an ordinary quark distribution function and a form 
factor 

H''{x,^,A') = d{x)F^{A^). (7) 

Here u{x) and d{x) are the unpolarized quark distributions and F'^^'^\A^) are defined through 
the electro-magnetic form factors of proton and neutron: 

= 2Ff + , F^ = Ff + 2Fi". (8) 

In the same context, the function H is written to be related to quark helicity distributions and 
axial- vector form factors: 

£"(a;,e,A2) = Auy{x)g\{A^)/gl{Q), 

H\x,^,A') = Ady{x)g'i{A')/gi{0), (9) 
where g'X = \gA + \g% gi = -|fi'A + and g^^ = ^gA- 

The functions E and E have no definite limit at ^ 0, ^ ^ 0, as it exists for the functions 
H and H (cf. Eq. (|l])). In absence of any other guide the ansatz for E is chosen in a form 
analogous to the function H: 

i?"(x,e,A2) = u{x)F^{A')/2, 

E\x,i,A^) = d{x)F,\A'), (10) 

where F^*-'^-* is defined in the same way as F^^'^^ in Eq. (§). 
The function E is modelled to be due to the pion pole: 

E^{x, e, A^) = -E^x, e, A^) = 1 E._,„,,(x, e, a^) , (n) 

E^.poieix, e, A^) = e{-^<x<0 hA{A') ^"^(fj ' (12) 
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where ^{z) = 3/4(1 — z'^) is the pion distribution amphtude, /i^(A^) = ^_^^2 , and 6{x) is the 
usual step function. 

To introduce, as a next step, a dependence of the GPDs on the skewedness parameter ^, the 
double-distribution formalism ^ can be used. In this model, the A^-independent part of the 
function H can be written in the following form 



Hhoi^, = f dy f ' dt6{x~y~ h{y, t) q{y) . (13) 
Here q{y) is the ordinary quark distribution and h{y, t) is the so-called profile function: 



22»+ir2(6+l) (1 - Ij/D^'+i ' 

with 6 as a free parameter and 6 — oo corresponding to the skewedness-independent parame- 
terization. Analogous expressions can be written for the functions H and E. 

The functions H and E in the form of double distributions lack the correct polynomiality 
properties of GPDs. However, they can be restored by introducing the so-called D-term ||14|| . 
The D-term contribution has different sign for H and E and hence cancels in Ji's angular 
momentum sum rule Eq. (^). The full model expressions for the GPDs therefore have the 
following form: 

= Kd{x,0 - Oi^-\x\)^D{^), (15) 

where Nj = 3 is the number of active flavors. The parameterization of the D-term is taken 
in a form that follows from chiral quark soliton model calculations |T^. For more details on 
possible GPD parameterizations and their discussion we refer to Ref.s |12, |13[. 

Quark number density and quark helicity distributions q{x) and Ag(x), respectively, were used 
in the parameterizations of Ref. ^ (at Q'^ = 2 GeV^) and Ref. (at = i QeV^). 

In the projections described below five different versions of GPD parameterizations were in- 
cluded: 

(A) H{x,^), H{x,^), and E{x,^) are skewedness independent and given by Eq.s (0), (|^), and 
([To|), respectively. 

(B) The skewedness dependence of H{x,^), H{x,^), and E{x,^) is generated through the 
double-distribution formalism Eq.s ([T3|) and ( p^lf ) with the parameter 6=1. 

(C) same as (B), but with the parameter b = 3. 

(D) same as (B), but the contribution of the D-term to H{x,C,) and E{x,C,) is included addi- 
tionally, according to Eq. (p!5|) . 

(E) same as (D), but with the parameter 6 = 3. 



For E(x,^) the pion pole expression Eq. ([Tl| ) is used for all versions of GPDs. 
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Figure 1: Illustration of the x-behaviour of various generalized parton distributions for different 
values of the skewedness parameter ^ (at A"^ = 0). For explanations see text. 

The dependences of the GPDs on the variables x and ^ (at = 0), obtained from these five 
different parameterizations, are displayed for a few examples in fig. Panels a) and b) show 
the GPD H in version (B) for different values of ^, compared between u- and d-quark. Panels 
c) and d) show the unpolarized u-quark GPDs H and E at fixed C, = 0.4, compared between 
different versions: (B) as solid line, (C) as dashed line and (D) as dotted line. Panels e) and 
f) show the polarized w-quark GPDs H and E in version (B) compared for different values 
of ^. Although E does not appear in the DVCS asymmetries (see below), its u- vs. (i-quark 
comparison is shown, as well. By definition E is large at = (see Eq. ([Tl|)). 



5 



3 Deeply Virtual Compton Scattering 



In the Bjorken limit, DVCS is dominated by the liandbag diagram (fig. 0) and its amphtude 
can be factorized into a soft part described by GPDs and a hard part representing a parton 
process calculable in perturbative QCD. In this limit the skewedness parameter ^ can be related 
to the Bjorken variable xb'- 

xb/2 



(16) 



1 - xb/2 

The same final state, epj, can also be produced via the Bethe-Heitler (BH) process in which an 
electron scatters elastically off the target proton and the initial or final state electron radiates 
a real photon. The cross section of this process can be calculated exactly once the Dirac 
and Pauli nucleon form-factors are known. On the one hand, the BH process constitutes the 
main background to DVCS, on the other hand their interference opens the unique opportunity 
for independent measurements of the real and imaginary parts of a certain DVCS amplitude 
combination (see below). 

The cross-section of the DVCS process and its interference with the BH process has been 
considered in a number of papers P, |^ . The kinematic configuration of the process ep — > epj is 
shown in fig. ^. Here 0^ describes the azimuthal orientation of the production plane (comprising 
7*, 7 and p) relative to the scattering plane (comprising initial and final lepton, as well as the 
virtual photon). The laboratory polar angle between virtual and real photon is denoted by 6.y^*. 
In-plane (i.e. 0^ = 0) differential cross-sections for DVCS, BH and total 7 production in e^p 
interactions at the HERMES energy of = 27.5 GeV were calculated following Ref. and 
are presented in fig. ^. The DVCS cross-section shows a maximum at G-y^* = 0, while the BH 
cross-section has a three-pole structure due to the propagators of the virtual electrons and the 
virtual photon. It is apparent that at HERMES energies the cross-section of DVCS is smaller 
than that of BH over practically the entire kinematic region. Therefore measurements of GPDs 
at HERMES must be based upon the interference between the DVCS and BH processes. 





Figure 2: Leading handbag dia Figure 3: Kinematic configuration for the process 
gram for DVCS. gp > 



In Ref. amplitudes of DVCS, BH and of the interference terms were calculated at the 



leading twist-2 level for polarized and unpolarized initial particles. This approach was used to 
determine the below presented projections for future measurements of DVCS-BH interference 
effects at HERMES. 
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Figure 4: Illustration of the behaviour of the differential in-plane cross-section as a function 
of the polar angle between the virtual and the real photon for DVCS (dashed lines), Bethe- 
Heitler (dotted lines) and total 7 production (solid lines) in e~^p interactions at HERMES 
energy — 27.5 GeV. Different panels are for different values of xb and . 

The primary interest of HERMES in the context of the planned recoil detector upgrade lies 
with measurements using an unpolarized proton target. In this case two different types of 
experiments are possible that will be giving insight into GPDs: 

I) Measurements of the lepton charge asymmetry with unpolarized leptons of either charge: 



dAa,h = da{e+p) - da{e-p) ~ cos(</)^) x Re ^F.Hi + J^i^i + ^2):^! - ^^2^i| (17) 

This asymmetry allows access to the real parts of the DVCS amplitudes TYi, TYi, and £1. 
II) Measurements of the beam helicity asymmetry using a polarized positron beam: 

d/^ULu = da{e^p) - da{e^p) ~ sin(0^) x Imj FiT^i + ^^^(^i + ^2)7?i - ^^2^i| (18) 



This asymmetry allows access to the imaginary parts of the same amplitudes. 

The imaginary and real parts of the DVCS amplitudes Hi and Hi are related to the GPDs H 
and H as follows: 

ImHi = -nJ2eliHi^,^,A')-Hi-^,^,A')), 
q 

ImHi = -nY.eliH{^,^,A^) + H{-^,^,A^)), 
q 

ReHi = j:el[pJ^'Hix,^,A'){-^ + -L^dx 



ReH 



where P denotes Cauchy's principal value. The DVCS amplitudes Si and Si can be expressed 
through E and E analogously. 

Projections for statistical accuracies attainable in measurements of the lepton charge and lepton 
beam helicity asymmetry at HERMES were calculated for an integrated luminosity of 2 fb~^ 0] 
which corresponds to the expected value for one year of data taking. The HERMES geometrical 
acceptance for the detection of the scattered electron, the photon and the recoil proton was 
taken into account. The following kinematic cutsQwere applied: Ee > 3.5 GeV, W"^ > 4 GeV^, 
g2 > 1 GeV^, E.^ > 1 GeV, Pp > 0.2 GeV, 0.35 < 6^^^ < 1.35 rad, and 15 < 9^^. < 70 mrad. 

To calculate the asymmetries Eq.s ([T7|) and (|18]) the 5-fold differential cross-section has to 
be integrated over the appropriate kinematic region accounting for the HERMES geometrical 
acceptance: 

dxdQ'^d\A^\d(f)^d(j)ei , 



u dxdQ'^d\A'^\d(j)^d(f)ei 

where (j)ei is the azimuthal angle of the scattered electron. It is appropriate to define differences 
and sums of certain cross-sections: 

dAach 



dEach 



dAa 



LU 



dT,aLu 



da{e^p) 


da{e p) 


d(f)^ 


d(f)^ 


da{e~^p) 


1 da{e-p) 


d(f)j 


d(j)^ 


da{e^p) 


da{e^p) 


d(j)^ 


d(j)^ ' 


da{e^p) 


1 da{e^p) 


d(j)y 


d(f)^ 



(20) 



Using these definitions the 0-dependence of the the lepton charge asymmetry reads 



AM) = tZZ^S^^'Z . (21) 



^Here and E.y are the energy of the incoming electron and the outgoing photon, respectively, while Pp is 
the momentum of the outgoing proton. The standard DIS variables and describe the four-momentum 
transfer and the energy of the 7*p system, respectively. The laboratory polar angle of the outgoing proton is 
given by 0'"''. 
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while an integrated lepton charge asymmetry can be defined by forming the difference between 
two integrals over appropriately defined 'halves' of the 0-distribution: 



^ _ d<P dAach/ d(j) - J^jj^ dcj) dAach/ dcj) ^^^^ 
Jq"" d4>dT,ach/d(f) 

Analogously the 0-dependent and an integrated lepton beam helicity asymmetry, respectively, 
are defined as follows: 



J^_Jd(j)dAaLu/d(l) 
U-A^ d(j)dJ:(TLu/d(l) 



^ £ dcj) dAam/d^ - dAam/d^ ^^4) 

J^'^ d(j) dJ^aiu I d(f) 

Projections for the statistical accuracy attainable in measuring the lepton charge asymmetry 
are presented in fig. ^ The results are shown as a function of the azimuthal angle 0^ in 
two regions of x^, using different GPD parameterizations. The asymmetry clearly depends on 
the particular parameterization and can even change its sign in dependence on Xb- Therefore, 
in the experiment the asymmetry will have to be studied differentially as much as possible. 
Apparently, the inclusion of the D-term leads to essential changes in the asymmetry and future 
measurements at HERMES may allow to confirm its importance experimentally. In fig. |^ the 
integrated lepton charge asymmetry Act is shown as a function of xb and A^, based on the 
GPD parameterization (E). No A^-dependence is seen on the basis of the chosen GPD model. 
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Figure 5: Projected statistical accuracy for the lepton charge asymmetry Eq. (|2T|) as a function 
of the azimuthal angle cj)^ between scattering plane and production plane. Predictions of different 
GPDs models (see text) are shown in two regions of xb- Version (A) - dash-dotted line; version 
(B) - long-space dotted line; version (C) - long-space dashed line; version (D) - dotted line; 
version (E) - dashed line. 



Projections for the statistical accuracy attainable in measurements of the helicity asymmetry 
are presented in fig. |^, as a function of the azimuthal angle (f)^ (left panel) and as a function of 
Xb and A^ (right panel). As can be seen, the projections of the statistical accuracy promise a 
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Figure 6: Projected statistical accuracy for the (p-integrated lepton charge asymmetry Eq. 
as a function of xb and based upon GPD model version (E). 



considerable improvement compared to the present HERMES measurement shown additionally 
in fig. 1^ (left panel). The projected asymmetry changes only slightly in dependence on the 
GPD parameterizations when the parameter h of the profile function Eq. (0) varies in the 
range (1 ; oo). Note that the beam helicity asymmetry is not sensitive to the D-term as it does 
not contribute to the imaginary parts of DVCS amplitudes. 




I 5 1 — , — , — , — , — I — , — . — , — , — I — , — , — , — , — I 
0.1 0.2 0.3 

Xb 

Figure 7: Projected statistical accuracy for the beam helicity asymmetry Eq.s 0) as a 
function of the azimuthal angle 0^ between scattering plane and production plane (left panel, 
closed circles) and as a function of Xb and (right panel). The line conventions are as in 
fig^. The present measurement of the azimuthal asymmetry by HERMES is shown in the left 
panel by open circles. 
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The considerations presented above are based upon the leading twist-2 level using the am- 
phtudes calculated in Ref. |18|. More elaborate analyses include twist-3 effects |jl5|, |19| and 



next-to-leading order calculations . Nevertheless, in the present situation where the general- 
ized parton distributions are practically unknown, the approach adopted in this paper appears 
adequate to evaluate the expected size of the predictions. 

It appears worth noting that additional important constraints on GPDs can be expected from 
the analysis of the same data set by studying hard exclusive production of both pseudo-scalar 
and vector mesons. As compared to DVCS these processes will provide information on different 
combinations of generalized parton distributions. 



4 Conclusions 



Expected statistical accuracies have been evaluated for future measurements of Deeply Virtual 
Compton Scattering with the HERMES spectrometer complemented by a to-be-built recoil 
detector. Using polarized electrons and positrons of HERA with different helicities, in con- 
junction with an unpolarized proton target, it becomes possible to measure the lepton charge 
asymmetry and the lepton beam helicity asymmetry, both induced by the interference of the 
DVCS and the Bethe-Heitler process. From these asymmetries the real and the imaginary part 
of a certain DVCS amplitude combination can be determined. 

The expected size of the asymmetries has been evaluated using various parameterizations of 
the underlying generalized parton distributions. A number of different parameterizations has 
been used to compensate as much as possible for the present poor knowledge on GPDs. The 
level of the attainable statistical accuracy is mainly determined by the cross section of the 
Bethe-Heitler process that dominates the reaction ep —>■ epy at HERMES energies. 

It has been shown that the planned measurements of hard exclusive real photon production at 
HERMES will be of high statistical significance. The envisaged separate results on the real and 
the imaginary part of a certain DVCS amplitude combination will constitute an important step 
towards the determination of the generalized parton distributions. The measured constraints 
will serve as very useful input for the further modelling of generalized parton distributions. 

We thank R. Kaiser for careful reading of the manuscript. 
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